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Abstract 

To discuss the interplay of electronic and lattice degrees of freedom in sys- 
tems with strong Coulomb correlations we have performed an extensive nu- 
merical study of the two-dimensional Holstein t-J model. The model de- 
scribes the interaction of holes, doped in a quantum antiferromagnet, with a 
dispersionsless optical phonon mode. We apply finite-lattice Lanczos diago- 
nalization, combined with a well-controlled phonon Hilbert space truncation, 
to the Hamiltonian. The focus is on the dynamical properties. In particular 
we have evaluated the single-particle spectral function and the optical con- 
ductivity for characteristic hole-phonon couplings, spin exchange interactions 
and phonon frequencies. The results are used to analyze the formation of 
lattice hole polarons in great detail. Links with experiments on layered per- 
ovskites are made. Supplementary we compare the Chebyshev recursion and 
maximum entropy algorithms, used for calculating spectral functions, with 

standard Lanczos methods. 
PACS number (s): 71.38.+i 



1 



Typeset using REVTgX 



I. INTRODUCTION 



There has been renewed interest in the analysis of strongly coupled elect ron-phonon (EP) 
systems since it was recognized that intrinsic polaron-like lattice distortions are a common 
feature of several important classes of perovskites such as the high-T c cuprates the 
non-metallic nickelates ||, and the colossal magnetoresistance manganites [|]||. In these 
compounds, e.g., in La2_ x Sr x [Cu/Ni]04 or Lai_ x Ca x Mn03, the charge carriers susceptible 
to lattice polaron formation mediate the antiferromagnetic (AFM) or ferromagnetic (FM) 
interactions between the Cu 2+ , Ni 2+ (S = 1/2, 1) or Mn 3+ (S = 3/2) ions (core spins), 
suggesting that the local lattice or Jahn-Teller distortions play also a significant role in 
determining the electronic and magnetic properties ||,|7]]. 

A number of experimental observations indicate spatial modulations of the spin and 
charge density in the copper-, nickel- and manganese-based perovskites. For example, 
dynamical stripe correlations of spins and holes were found from neutron scattering || and X- 
ray absorption (EXFAS) || in superconducting Lai.gsSro.isCuCU and YBa2Cu3C>6.6- Rather 



static polaronic superlattice structures have been detected by electron [10], neutron [ I I 



and X-ray diffraction QT2j| measurements in hole doped (insulating) La2_ x Sr x Ni04. Local 
structure studies of Lai_ x Sr x Mn0 3 , performed using pulsed neutron powder diffraction and 
atomic pair-density function analysis []I3"| |, suggest that nonuniform charge distributions 
related to lattice polarons exist in the paramagnetic, AFM and FM phases, and are even 
present in the metallic state up to x ~ 0.35 . 

A striking feature of all these charge and spin ordering phenomena is their cooperative 
nature, demonstrating the strong interaction between spin, charge and lattice degrees of 
freedom |6|JI^[. Accordingly, the lattice anomalies are well reflected in the transport prop- 
erties. For example, there is observed a striking suppression of superconductivity in the 
[La, Sr^CuC^ family compounds with ~ 1/8 doping, which is believed to have a connection 
with charge and spin stripe formation ||. Undoubtedly, structural correlations and (po- 
laronic) charge order are shown to be the driving force behind the pronounced resistivity 



anomalies found in [La, Sr^NiC^ for hole concentrations 1/3 and 1/2 fl0|]T5|j . 

Since polaronic effects might be responsible for some of the puzzling thermodynamic and 
transport properties of the cuprates, nickelates and manganites, increasing effort has been 
aimed to test the polaronic nature of charge carriers by means of optical experiments. One 
of the physical quantities which contains valuable information about the low-energy ex- 
citations in interacting EP systems is the optical conductivity, ct{uj), usually determined 
from reflectivity measurements. Here the anomalous broad optical absorption detected 
for the La2_ x Sr x [Cu, Ni]04 compounds in several (normal-state) infrared experiments has 
attracted a lot of attention. In particular for the nickelate system, where the existence 
of self-trapped charge carriers is well established, the mid-infrared absorption band has 
been interpreted successfully in terms of photon-assisted hopping of small-size lattice po- 
larons |16| , |17[ . On the other hand, from a comparison of the optical absorption of several per- 
ovskites (Lai.67Sr .33NiO4, Sri.sLao.sMnOzi, and Ndi.geCeo.o/iCuOzi), it has been argued that 



in the low-doping cuprates the lattice polarons appear to be rather large JT8|,|T9[ . Although 
recent time-resolved optical relaxation measurements on YBaCu07_,5 give additional evi- 
dence for the existence of polaronic carriers in the cuprates [|2U|] , it is worth mentioning that 
the present experimental situation certainly does not permit to draw any decisive conclusion 
concerning the relevance of the polaron/bipolaron scenario for high-T c superconductivity. 

From a theoretical point of view, the main problem for a proper treatment of such sys- 
tems, exhibiting besides a substantial EP coupling strong Coulomb correlations, is that the 
coherent motion of charge carriers, heavily dressed, e.g., by interaction with magnons, takes 
place on a strongly reduced energy scale being almost comparable to the relevant phonon 
frequencies. As a result the standard strong-coupling Migdal-Eliashberg approach based on 
the adiabatic Migdal theorem might break down PH, 22], and it has been argued that non- 



adiabatic vertex corrections are responsible for the enhancement of T c p3|-t25|. Furthermore, 
as a consequence of the pre-existing magnetic "self-localization" of the charge carriers, pola- 
ronic effects play an important role even at rather moderate EP coupling strengths 



Therefore the investigation of such highly correlated EP systems is of fundamental impor- 



tance, not only in connection with the high-T c problem. 

Our aim here is to study what is perhaps the minimal microscopic model with respect to 
a strong spin-charge-lattice coupling, the so-called Holstein t-J model (HtJM). The focus 
is on the spectral properties of the HtJM Hamiltonian in the single-hole subspace relevant 
to the physically most interesting low doping regime. In order to obtain reliable results for 
the complete range of exchange interactions, EP couplings and phonon frequencies we will 
employ numerical exact diagonalization (ED) techniques on finite lattices, using massive 
parallel computers such as the CRAY T3E. 

The paper is organized as follows: In the next section we introduce the HtJM and 
comment on previous numerical investigations. Sec. Ill presents our ED results for the 
single-particle spectral function, which will be used in Sec. IV to discuss the hole-polaron 
band formation. In Sec. V we analyze the optical response in the HtJM. Our main results are 
summarized in Sec. VI. A comparative study of several computational techniques proposed 
for calculating spectral properties of very large Hamiltonian matrices is performed in the 
Appendix. 

II. THE HOLSTEIN T-J MODEL 

Our model is the standard two-dimensional t-J model appended by an additional 
Holstein-like interaction term with a dispersionsless (optical) phonon branch: 

W = -*E (4^ + H.c.) + J£ (S t S 3 - -n^) 

(yV (ij) 

-y^^E (6t + h) K + ( 6 & + \) • (i) 

i i 

The Holstein t-J Hamiltonian (JJ) acts in a projected Hilbert space without doubly occupied 
sites, where e£: = c£?(l — h i _ a ) is a spin-a electron annihilation (creation) operator on 
Wannier site i, hi = X)<j $i<fiia-> an d Si = J2 a ,a' Ci a T aa ,c ia , . Within an effective single-band 
description, the parameters t and J measure the transfer amplitude and the antiferromag- 
netic exchange interaction between nearest-neighbour sites on a square lattice, where t > J 
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corresponds to the situation in the cuprates. For example, J/t ~ 0.4 with t ~ 0.3 eV is 
commonly used to model the La2_ x Sr x Cu04 system. In ([[]), we have chosen the coupling 
to the hole, hi — 1 — hi is the local density operator of the spinless hole, as the dominant 
source of the EP coupling, because in the t-J model the unoccupied site corresponds to the 
Zhang-Rice singlet (formed from Cu3d x 2_ y 2 and 2p xy hole orbits) for which the coupling 
should be much stronger than for the occupied (Cu 2+ ) site [|3~l"H3~3"l ■ The hole-phonon cou- 
pling constant is denoted by e p , and hu>o is the bare phonon frequency (below h — 1 and 
all energies are measured in units of t). Then g« = {b\ + b j )/^2MojQ may be thought of 
representing an internal vibrational degree of freedom of the lattice site i (bf^ annihilates 
(creates) an optical phonon). Indeed for La2[Cu, Ni]C>4 the oxygen vibrational modes have 
a small dispersion in the Brillouin zone, i.e., they are revealed to be very localized in real 
space and can be considered as independent [j34|j . 

Whereas the electronic properties of the pure t-J model have been studied extensively in 
the context of high-T c superconductivity JJS], very little is known theoretically about the in- 
terplay between magnetic exchange interaction and EP coupling in the framework of the 2D 
HtJM. The changes of the quasiparticle (QP) properties due to the cooperative/competitive 
effects of the hole-phonon and hole-magnon interactions are expected to be very complex 
and as yet there exist no well-controlled analytical techniques to solve the problem partic- 
ularly in the physically relevant coupling regime, where the characteristic electronic (t, J) 
and phononic (s p , u ) energy scales are not well separated. Naturally such a dressed hole 



QP will show the characteristics of both "lattice" and "magnetic" (hole) polarons |36 



Early attempts to study lattice polaron effects in the HtJM were based on a variational 



squeezed-polaron approach to the phonon subsystem ||30|| . The effective electronic Hamilto- 
nian, obtained in the transformed phonon vacuum state, was treated by means of Gutzwiller 
or slave boson techniques and the main features of the ground state phase diagram of the Hol- 



stein Hubbard/t-J models were worked out on a mean-field level of approximation |37| , |3"8 
Analyzing the static charge susceptibility of the HtJM, a transition to a Peierls distorted 
phase was found to occur at quarter filling ||39|| . In the low-doping regime and for moderate 



EP interactions the lowest polaron-magnon band states were determined for the Lang-Firsov 
transformed 2D HtJM within the zero-phonon and spin-wave approximations, indicating 
the formation of large-size (mainly magnetic) hole polarons with soliton-like-shaped highly 
asymmetric wave functions, oriented along the diagonals of the square lattice P0| . Com- 
bining the spin-wave approximation with an iterative Lanczos algorithm, the ground-state 
properties of the 2D HtJM were examined including multi-phonon states f41 ]. A sharp 



crossover to small-size less-dispersive symmetric lattice hole polarons was obtained at a 
critical value of the EP coupling fril l. Using the self-consistent Born approximation, the 



phonon-induced mass renormalization of a hole in the HtJM was found to depend strongly 



on the confining string potential due to the spin background [^T]]. Whereas a large AFM 
exchange interaction certainly enhances the tendency towards lattice polaron formation, 
there is evidence that in the limit J — > the influence of phonons becomes even smaller 



than in the corresponding uncorrelated (spinless fermion) model ||21|| . More recently, the 



single-hole spectral density function and the optical conductivity were studied within the 
non-crossing approximation for spin-wave and phonon interactions in a Holstein t-J-like 



Hamiltonian [22]. Perhaps surprisingly, the numerical estimate of the lowest order vertex 
corrections to the hole-phonon coupling vertex gives relatively small values, indicating that 
Migdal's approximation can be used in calculating the hole self-energy. This effect has 
been attributed to the significant suppression of the QP residue for a correlated electron 



system j22 



Of course these analytical approaches provide some insight into the problem, but it is 
difficult to judge how close to the actual properties of the HtJM the results are. From 
this point of view it is of interest to carry out exact diagonalizations of the HtJM on small 
clusters. Along this line, a number of ED studies have been performed treating the phonons 
in the static limit, i.e., classically, in the so-called adiabatic or frozen phonon approxima- 



tion |42|J3T| , |43|J^] . The stability of the ground state against static lattice distortions was 
investigated and the phase diagram of the 2D HtJM has been explored in the one- and 
two-hole sectors 1 31J1 , as well as for the quarter-filled band case [Q, where a polaronic 



superlattice was found to occur above a critical EP coupling strength. Neglecting, how- 
ever, the phonon dynamics, certainly overestimates the tendency towards lattice polaron or 



charge-density-wave (CDW) formation p5| , ^6| . As a first step to overcome the static treat 



ment of the Ht JM, a variational Lanczos diagonalization technique has been proposed |f26[ . 
which is based on an inhomogeneous modified Lang-Firsov transformation and allows for 
the description of static displacement field, non-adiabatic lattice polaron and squeezing ef- 
fects on fairly large clusters. Extending the ED approach to the full quantum phonon model 
requires a truncation of the bosonic part of the Hilbert space, which, strictly speaking, is 
infinite dimensional. Including at each site besides the phonon vacuum a second (coherent) 
displaced oscillator state, the ID HtJM has been analyzed with respect to the occurrence 
of CDW order, superconducting fluctuations, and phase separation [|45|| . Whether the com- 



plex (phonon-assisted) transport in the HtJM can be described with sufficient accuracy 
within such a simple two-state variational approach remains questionable, at least in the 
intermediate coupling regime. For the 2D case more advanced truncation procedures in 
momentum space, retaining a single or small number of phonon modes pi|,|?7| , have been 
employed to show that dynamical in-plane phonons could help to stabilize inhomogeneous 
stripe-like phases. Taking into account all dynamical phonon modes, in a previous work E7 



the authors have examined the ground-state properties of the 2D Holstein t-J model in the 
one- and two-hole subspaces by means of direct Lanczos diagonalization. Most notably, the 
transition from delocalized to self-trapped lattice hole polarons (hole bipolarons), signaled 
by a suppression of the kinetic energy and an enhancement of the hole-phonon (hole-hole) 
correlation functions, was found to be favored by strong Coulomb correlations. 

In the present paper, we extend this rigorous approach, which preserves the full dy- 
namics and quantum nature of phonons, to the calculation of dynamical quantities. In 
particular we compute the single-hole spectral function and the optical conductivity to dis- 
cuss the transport of hole polarons in the framework of the 2D HtJM. The technical details 
of the computational procedures such as the basis implementation and symmetrization, the 
phonon Hilbert space truncation and the spectral moment approach will be discussed in the 



Appendix. 



III. SINGLE PARTICLE SPECTRAL PROPERTIES 

In the numerical analysis of the 2D HtJM we start with a discussion of the single-hole 
spectral function 



n,cr 

Figure 1 displays A^(o;) for the allowed values of the momentum K of a ten-site square 
lattice. To visualize the intensities (spectral weights) connected with the various peaks 
(excitations) in each K-sector we also have shown the integrated density of states 

N(u)= f dJjjY, A *W- (3) 

°° K 

Of course, in the absence of EP coupling we reproduce the single-particle spectrum of the 
pure t-J model, which for comparison is shown in Fig. 1 (a). In this case the single-hole 
excitations are reasonably well understood |35|| . Ag(u) contains a "quasiparticle" peak 
at the bottom of the spectrum followed by a lump of spectral weight with some internal 
structure. The QP pole, corresponding directly to the coherent single-hole ground state 
with momentum (37r/5, 7r/5), is separated by a pseudogap of approximative size J from the 
lower edge of a broad incoherent continuum 6t wide. Despite the significant redistribution 
of intensity to higher energies as a consequence of multiple spin-wave (spinon) scattering 
of the magnetic polaron [48j, the QP residue still contains a substantial amount of spectral 



weight. It has been argued, that the QP peak and the low-energy satellite structures, 
which can be identified as the next excited states of the string problem, will survive in 
the bulk limit ||35|| . The momentum dependences of Ag(u) indicate that the low-energy 
structure of the i^-integrated spectrum is dominated by excitations with wave vectors near 
the magnetic Brillouin zone boundary (particularly by the i^-vectors closest to (7r/2,7r/2) 
and (it, 0)) [|49|| . On the other hand, for momenta K = (0, 0) and (tt, tt) most of the weight 
of is concentrated at high energies, i.e., in the incoherent part of the spectrum. 
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Now we consider the influence of an additional weak EP interaction (Fig. 1 (b)). In 
the weak-coupling regime the mass renormalization of the coherent QP band due to the 
hole-phonon coupling is generally small compared with that arising from the hole-spin 
interactions (magnetic polaron regime). In particular, the integrated density of states are 
barely changed from those of the pure t-J model. The new structures, nevertheless observed 
in the Ag spectra shown Fig. 1 (b), correspond to predominantly "phononic" side bands 
separated from the particle-spin excitations by multiples of the bare phonon frequency ujq. 
However, these phonon resonances have less and less "electronic" spectral weight the more 
phonons are involved. This is because Ag(u>) measures the overlap of these excited states 
with the state obtained by creating a hole in the zero-phonon Heisenberg ground state. 

To analyze the formation of lattice hole polarons in the 2D HtJM, we have compared 
the non-adiabatic weak- and strong-coupling situations in Fig. 2. Looking first at the case 
J = 0.4, panels (a) and (c), we observe that with increasing e p the lowest peaks in each 
Ag- start to separate from the rest of the spectrum. These states become very close in 
energy and finally a narrow well-separated lattice hole-polaron band evolves in the strong- 
coupling case (e p = 4.0). Since the gap to the next energy "band" is of the order of u (see 
panel (c)), these excitations will be triggered by an one-phonon absorption process. The 
transition to the lattice hole-polaron state is accompanied by a strong increase in the on-site 
hole-phonon correlations |27|], indicating that the lattice polaron QP comprising a "quasi- 
localized" hole and the phonon cloud is mainly confined to a single lattice site (small-size 
lattice hole polaron). Compared to the non-interacting single-electron (Holstein model) |5(J 
or spinless fermion (Holstein-t model) |26| problems, the critical EP coupling strength for 
lattice polaron formation is considerably reduced due to magnetic prelocalization effects. 
The phonon distribution function shows that the small-size lattice hole-polaron state is 
basically a multi-phonon state (see Ref. [[51]; cf. also Fig. 8 in the Appendix). For example, 
at e p = 4.0 and J = 0.4 we found the mean phonon number in the ground state to be about 
4-5. Here the phonons will heavily dress the hole and the QP pole strength becomes strongly 
suppressed (cf. the discussion in Sec. IV). At the same time spectral weight is transferred to 
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the high-energy part and the whole spectrum becomes incoherently broadened. Therefore, 
we observe an overall smoothing of N(u). We note that the vertex corrections to the 
hole-phonon vertex are non-negligible for the large hole-phonon couplings and intermediate 
phonon frequencies considered here (e p , y/e p uj ^> J). This has to be contrasted with the 
weak to moderate EP coupling regimes, where the vertex correction was found to be much 
smaller than unity [^2[ . 

Next, we want to discuss the effect of the exchange interaction J. For the t-J model it is 
well known that as long as J/t is large, i.e., the spins are antiferromagnetically ordered, the 
derealization of a hole is hindered as a consequence of a linear confining string potential 
originated from the misaligned spins along the hole hopping paths. If J tends to zero, 
the string potential is weakened and the hole performs large excursions polarizing the spin 
background ferromagnetically. In a finite system this so-called Nagaoka transition takes 
place at a finite critical exchange interaction (J^ = 0.13807 for the ten-site lattice). For 
J < Jn, the single-hole ground state has momentum K = (0, 0) and is fully spin aligned. 
At J = 0, the spectrum is completely incoherent, extending over almost the unperturbed 
bandwidth 8t. As can be seen from Fig. 2 (b), showing the case J = 0.1 < J^, a weak 
EP coupling does not affect this scenario. Quite the reverse, from a self-consistent Born 
calculation it was deduced that the influence of phonons on such a "Nagaoka spin-polaron" 
becomes even smaller than in the free fermion case f2l |. The situation, however, is changed 



completely in the strong EP coupling regime. As e p increases, the mobility of the hole is 
reduced and finally a strong local lattice distortion traps the hole almost on a single site. 
The hole-trapping makes even a small AFM exchange interaction very effective and we 
observe, as for the static (u — > 0) HtJM P5fl , an enhancement of the AFM correlations in 
the spin background. Consequently, the whole spectrum looks very similar to that obtained 
for J = 0.4 (cf. panels (c) and (d)). The only difference is that now the first spin excitations, 
having an excitation energy oc j( 2 / 3 ) ||52j| , occur just above the lattice hole-polaron band, 
thus filling the gap between the first and second polaronic bands. 

Finally, we consider the regime, where the bare phonon frequency is much larger than 
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the effective hole bandwidth of the t-J model. For 0.1 < J < 0.4, the coherent bandwidth 
can be approximated by AE l ~ J ~ 1.5J a with a ~ 1 Figure 3 shows the results for the 



spectral function Ag{ui) at u; = 3.0, £ p = 4.0, and J = 0.4. In this case the lowest (QP) 
excitation band is only weakly renormalized by the phonons which can follow the motion 
of the magnetic polaron instantaneously. The weights of the first m-phonon states in the 
ground state are 0.436, 0.332, 0.15, and 0.055 for m = 0, 1, 2, and 3, respectively. As 
can be seen by comparing Figs. 1 (a), 2 (c) and 3, the QP pole strength is determined by 
the subtle interplay of e p and ujq. In the anti-adiabatic strong-coupling limit, the high- 
energy part of the Ag(uj) spectrum is clearly dominated by interband transitions according 
to multi-phonon absorptions (cf. the pronounced jumps in N(u) at energies w~nx u> ). 

IV. POLARON BAND FORMATION 

To illustrate the formation of the lattice hole-polaron band in some more detail, we have 
calculated the "coherent" band dispersion of the 2D HtJM, Eg, on a 16-site lattice. Eg was 
determined from the lowest pole of each Ag(u) having finite spectral weight. The resulting 
band structure is shown in Fig. 4 along the principal directions in the Brillouin zone. For 
the pure t-J model, the dispersion of the QP band is highly consistent with a simple tight- 



binding form |54| that includes hole hopping processes to first and second-nearest neighbours 
on the same sublattice only (leaving the magnetic order undistorted) . From Fig. 4 (a) the 
minima of the QP dispersion are found to be located at the momenta K = (±7r/2, ±7r/2) (the 
hidden symmetry of the 4x4 cluster leads to an accidental degeneracy with the K = (±7r, 0), 
(0, ±7r) states). Obviously, the energy dispersion is not significantly changed including a 
weak EP coupling, provided that the phonon frequency exceeds the effective bandwidth of 
the magnetic polaron (ujq = 0.8 > AE l ~ J ). On the other hand, in the adiabatic regime 
(c^o <§C AE l ~ J ), a remarkable "flattening" of the band structure of the HtJM takes place, 
which can be attributed to the "hybridization" of the hole QP band with the dispersionsless 
optical phonon. As a result the coherent bandwidth AE = supgEg —migEg is roughly 
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given by ojq. 

In order to make this discussion more quantitative, we have evaluated numerically the 
spectral weight of the different band states, 



X 0,K 1 Q-K,a< 0,Q I 



\ 0,Q ' Q—K,a Q—K,a> 0,Q ' 

where I^q^™ 1 '*) denotes the one-hole (polaron) state being lowest in energy within the Q- 
sector. The results for Zg are displayed in Fig. 4 (b). As already discussed in Sec. Ill, 
Zg is strongly reduced from unity in the pure t-J model (but seems to remain finite in 
the thermodynamic limit if J > 0). Including a weak adiabatic EP interaction (e p = 0.2, 
Uq = 0.2), the band states with energies larger than E + ujq change their character from 
"hole-like" to "phonon-like" . Consequently, we observe a nearly complete suppression of 
Z(q ,0) and Z^^/2) ■ Increasing the phonon frequency, this tendency is reversed (see the curve 
for ujq = 0.8). Note that similar effects have been discussed quite recently in terms of the 



2 ■ 



(4) 



single-electron Holstein model 



To complete this section, we show in Fig. 5 the variation of Zg with the EP coupling 
strength for the i?-vectors of the ten-site lattice studied in Sec. III. The wave-function 
renormalization factor Zg can be taken as a measure of the "contribution" of the hole 
(dressed at e p = by spin-wave excitations) to the polaronic QP (having total momentum 
K). The data obtained at weak EP coupling unambiguously confirm the different nature 
of band states in this regime: we found practically zero-phonon "hole" states at the band 
minima (K = (37r/5, 7r/5), triangles down) and "phonon" states, which are only weakly 
affected by the hole, around the (flat) band maxima (K = (vr,7r), triangles up). With 
increasing e p , a strong "mixing" of holes and phonons takes place, whereby both quantum 
objects completely loose their own identity. Concomitantly Z^ decreases for the "hole-like" 
states but increases (first of all) for the "phonon-like" states. At large e p , a small lattice 
hole polaron is formed, which, according to the numerical results, has an extremely small 
spectral weight. Hence the question arises whether the lattice hole polaron is a "good" QP 
in the sense that one can construct a QP operator, — > dg a , having large spectral weight 
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at the lowest pole in the spectrum. Indeed, it was demonstrated in recent ED work that 
it is possible to construct such a composite electron/hole-phonon (polaron) operator by an 
appropriate phonon dressing of c Ka for the Holstein model |5(| as well as for the t-J model 



coupled to buckling/breathing modes p9 



V. OPTICAL CONDUCTIVITY 



Using standard linear response theory, the Kubo formula for the real (absorbent) part 
of the frequency-dependent optical conductivity gives two physically distinct contribu- 
tions 



u > 0. 



(5) 



The first so-called Drude term at uj = is due to the free acceleration of the charge carriers 
by the electric field and the second term, frequently called the "regular term" , is due to finite 
frequency dissipative optical transitions to excited quasiparticle states. More explicitly, the 



regular part can be written in spectral representation at T = as |5H] 



<C(") 



N + n 



E 



(N-l) 



E 



(N-l) 
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u-(E 



(N-l) 
n.K 



E 



0,K ' 



(6) 



n,K 0,K 

where the summation is taken over the complete set of eigenstates with excitation energies 



[E 



(N-l) 
n.K 



E 



(N-l), 
0,K J 



in the one-hole (N — 1 electron) subspace. For the HtJM the 



(paramagnetic) current density operator has the form 



(7) 



Actually in (|]) an optical transition can take place only within the A'-sector of the ground 
state. 

The expression (||) has been used to calculate a re9 {uS) numerically for the t-J (e.g., see 
Ref. fl35| ) and Holstein models |59|j60|j50[1 on finite lattices employing ED techniques. First 



results, obtained recently for the ID HtJM [32], indicate that EP coupling effects may be at 



13 



least partly responsible for the experimentally observed frequency dependence of the optical 
conductivity in the CUO3 chains of the 1:2:3 family of high-T c oxides. 

In Fig. 6 we show a re9 (u) of the 2D ten-site HtJM, determined from (0) with at most 
M — 15 phonons, for typical values of the hole-phonon coupling and J = 0.4. In the weak 
EP coupling regime and for phonon frequencies ujq ~ AE l ~ J , we recover the main features 
of the optical absorption spectrum of the 2D t-J model i.e., an "anomalous" broad 



mid-infrared band [J ~ uj ~ 2t], separated from the Drude peak [D5(u); not shown] by a 
"pseudo-gap" ~ J, and an "incoherent" tail up to uj ~ It (compare panels (a) and (b)). 
Note particularly the even quantitative agreement between the a;-integrated spectral weight 
functions 



S"»= / du'aZ 9 ^') (8) 
Jo 

depicted in Figs. 6 (a) and (b). Since for the parameters used in (b), e p = 0.1 and ujq = 
0.8, the averaged kinetic energy per site (—7i t /N) is only weakly renormalized by the EP 
interaction (cf. Fig. 10 of Ref. [Z7||; 7i t denotes the first term in ((]])), also the Drude weight 



ire 2 



V = —(-H t )-2S^(oo), (9) 

obtained via the f-sum rule, should approximately be the same as for the non-interacting 
system (s p = 0). As discussed in the preceding section, the coherent band structure gets 
stronger renormalized at large e p , i.e., the QP band states have less spectral weight. At 
the same time the phonon distribution function in the ground state, |c™| 2 (see Eq. (18) 
Appendix), becomes considerably broadened p7|j51| , and the position of its maximum is 
shifted to larger values of m. Consequently the overlap with excited multi-phonon states 
is enlarged and the optical response is enhanced at higher energies. This redistribution of 
spectral weight from low to high energies can be seen in Fig. 6 (c). As expected the transition 
to the lattice hole-polaron state, at about e c (J = 0.4, ujq = 0.8) ~ 2.0, is accompanied by 
the development of a broad maximum in a re9 (u), whereas the Drude weight as well as 
the low-frequency optical response become strongly suppressed. Contrary, in the strong- 
coupling anti-adiabatic regime (e p = 4.0, uj = 10.0; see panel (d)), the low-frequency part 
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of a reg (u) is much less affected, but the optical spectrum shows additional superstructures 
corresponding to "interband" transitions between t— J— like absorption bands with different 



numbers of phonons |62] 



Let us now make contact with the experimentally observed characteristics of the mid- 
infrared spectra in the doped perovskites (for an review of the experimental and theoretical 
work on the optical conductivity see Refs. |63| , |35|j59| ,3| ) . Needless to say, that it is out of the 
scope of our ten-site cluster diagonalization study to give a quantitative theoretical descrip- 
tion of the complex optical properties of particular copper and nickel oxides. But we would 
like to stress that the 2D HtJM seems to contain the key ingredients to reproduce, at least 
qualitatively, the principal features of the optical absorption spectra of these compounds. 
This can be seen by comparing the results shown in Figs. 6 (b) and (c), which correspond 
to the weak and strong EP coupling situations realized in the cuprate (La 2 _xSr x Cu0 4 ) and 
nickelate (La2- x Sr x Ni04) systems, respectively. The "effective" EP interaction e p /t in the 
nickelates is estimated to be about one order of magnitude larger than in the cuprates simply 
because of the much smaller transfer amplitude (t ~ 0.08 eV $M). According to the inter- 



nal structure of the low-spin state, the hopping transport of spin-1/2 composite holes in a 
spin-1 background is rather complex; implying, within an effective single-band description, 



a strong reduction of the transition matrix elements |I4| . A striking feature of the absorption 
spectra in the cuprate superconductors is the presence of a broad mid-infrared (MIR) band, 
centered at about 0.5 eV in lightly doped La 2 _ x Sr x Cu0 4 (which, using t ~ 0.3 eV, means 
that ui ~ 1.5). Such a strong MIR absorption is clearly observed in Fig. 6 (b), which refers 
to the weak EP coupling case. Since it also appears in the pure t-J model (cf. panel (a)), 
this MIR band seems to be caused by the spin fluctuations around the charge carrier p5| . 
Obviously, it is quite difficult to distinguish the spectral weight, produced by the dressing 



of the hole due to the "bag" of reduced antiferromagnetism in its neighbourhood |S5 |, from 
other (e.g. hole-phonon coupling) processes that may contribute to the MIR band observed 
experimentally. The results presented for the HtJM in Fig. 6 (b) support the claims, how- 
ever, that the MIR band in the cuprates has a mainly "electronic" origin, i.e., the lattice 
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polaron effects are rather weak. The opposite is true for their isostructural counterpart, the 
nickelate system, where the MIR absorption band has been ascribed by many investigators 



to "polaronic" origin fllbjr8| . Within the HtJM such a situation can be modeled by the 
parameter set used in Fig. 6 (c). If we fix the energy scale by t ~ 0.08 eV (which is the 
estimate for the nearest neighbour transfer integral in La2_ x Sr x Ni04 given in Ref. [HJ), the 
maximum in the optical absorption is again located at about 0.5 eV. The whole spectrum 
clearly shows lattice polaron characteristics, where it seems that the lattice hole polarons 
are of small-to-intermediate size ||66|| . Most notably, we are able to reproduce the experi- 
mentally observed asymmetry in the shape of the spectrum, in particular the very gradual 
decay of a re9 {uo) at high energies. It is worth mentioning that this behaviour cannot be 
obtained from a simple fit to the analytical expressions derived for the small polaron hop- 
ping conductivity |67|j68|JT6[] . Exploiting the f-sum rule we found that there are almost no 



contributions from band-like carriers in agreement with the experimental findings flT6|,|18 . 

Finally, we would like to consider the case J = 0.1, which certainly has no relevance with 
respect to the perovskite materials but this parameter regime is interesting from a theoret- 
ical point of view. For e p = (t-J model), the momentum of the Nagaoka ground state 
is K = (0, 0) and the hole can propagate without any spin scattering in the ferromagnetic 
background. Therefore all the spectral weight stays in the Drude pole at uo = and the 
regular (incoherent) part of the optical conductivity vanishes exactly. In principle, at arbi- 
trary small EP interactions an optical transition can be achieved by "adding" phonons with 
momentum (-K) to band states Eg (in order to reach the K = (0,0) sector). However, 
the overlap of these excited states with the almost zero-phonon ground state is extremely 
small and, as a direct consequence, S re9 (oo) is negligible (see Fig. 7 (a)). A completely 
different situation arises when the charge carrier becomes "self-trapped" due to the lattice 
distortions. When the lattice hole-polaron formation takes place, the momentum of the 
ground state of the 2D ten-site HtJM is switched from (0,0) to K — (37r/5, 7r/5) and, what 
is more important, the AFM exchange interaction becomes more effective. This gives rise 
to the pronounced low-frequency peak structures found in Fig. 7 (b), which are originated 
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from spin excitations on an energy scale of the order of J. At higher energies we recover the 
characteristic features of the lattice hole-polaron spectrum discussed in Fig. 6 (c). 

VI. SUMMARY 

We have investigated the spectral properties of the two-dimensional Holstein t-J model 
as a generic model for studying polaronic effects in systems with strong Coulomb correlations. 
The use of purely numerical techniques allows to treat the electron and phonon degrees of 
freedom on an equal footing. In the numerical work the full dynamics of all the phonon modes 
was taken into account on finite lattices with up to 16 sites. The efficiency and accuracy of the 
employed phonon Hilbert space truncation, Lanczos diagonalization, Chebyshev recursion 
and maximum entropy methods have been demonstrated in the Appendix. 

In conclusion, we have shown that an additional hole-phonon coupling strongly affects 
the nature of the quasiparticle excitations and the charge-carrier transport in the 2D t-J 
model. Our results indicate that the magnetic prelocalization of the holes, due to strong 
spin correlations, strengthens the effect of the hole-phonon interaction. Thus the tendency 
towards lattice (hole) polaron formation is enhanced in strongly correlated electron systems. 

The wave-vector resolved single hole spectral function yields a flattening of the coherent 
band dispersion in the vicinity of the band maxima for the adiabatic weak and intermedi- 
ate hole-phonon coupling cases. As the hole-phonon interaction increases, a narrow lattice 
hole-polaron band separates from the rest of the spectrum but the spectral weight of these 
quasiparticle band states becomes more and more suppressed. In the anti-adiabatic regime 
the low-energy excitations of the t-J model are less affected by the hole-phonon interac- 
tion; at higher energies we found "interband" transitions corresponding to multi-phonon 
absorption processes. 

The formation of lattice hole polarons is accompanied by significant changes in the 
optical response of the system. Both the coherent Drude part and the low-frequency regular 
contribution to the optical conductivity are strongly reduced and, concomitantly, we observe 
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a substantial spectral weight transfer to higher energies. In the intermediate-to-strong 
hole-phonon coupling regimes, i.e. in the vicinity of the lattice hole-polaron transition, the 
lineshape of the optical absorption spectra is highly asymmetric with a long incoherent tail 
at high frequencies, which also seems to be in accordance with recent mid- and far-infrared 
optical experiments on La 2 _ x Sr x Ni04 pHyiSfl . The present calculations point toward the 
importance of strong spin correlation and polaronic effects in explaining the origin of the 
anomalous mid-infrared absorption band observed in the cuprate and nickelate compounds, 
respectively. Of course, for a realistic modelling of the complex high-spin low-spin nickelate 
system, one has to take into account more complicated hopping and spin exchange processes. 
Work along this line is in progress. 
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APPENDIX: NUMERICAL METHODS 

Symmetrized basis and Hilbert space truncation 

1. Basis symmetrization. The total Hilbert space of the HtJM can be written as the 
tensorial product space of electrons and phonons, spanned by the complete basis set 

{\<Puv} = \u)el®\v) ph } (10) 

with 
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N 

\ u )ei = n n (4r^ iou n iaiU e {o, i> , (n) 

TV 1 

l^ = n^^( & i) mi " l°U mi,«e{0,...,oo}. (12) 

i=l yJmi )V \ 

Here it = 1, . . . , D e \ and v — 1, . . . , _D p h label the basic states of the electronic and phononic 
subspaces with dimensions D e \ = (^)( JN ^~ JV<T ) and D ph = oo, respectively. 

Since the Hamiltonian (1) commutes with the electron number operator M e i = X^i(^i,t + 
h i: i) and the ^-component of the total spin S z = | Y^iLxi^i,] ~ the basis (pi]) has been 
constructed for fixed N h = N — N el and S z = S^ in . To further reduce the dimension of the 
total Hilbert space, we can exploit the space group symmetries [translations (G*t) and point 
group operations (Gx)] and the spin-flip invariance [(Gs); S z = O-subspace only]. Clearly 
working on finite bipartite clusters (N = k 2 + I 2 , k and I are both even or odd integers) 
with periodic boundary conditions, we do not have all the symmetry properties of the 
2D square lattice. Restricting ourselves to the one-dimensional non-equivalent irreducible 
representations of the group G(K) = Gt x Gl(K) x Gs, we can use the projection operator 

^,„ = -^ £ xf„e (is) 

yv-'V g^G(K) 



(with [H,Vg jrs ] = 0, Pt rs = V Rrs and 7^ irs 7^,^, = Vg rs Sjtf, 5 r y 5 S)S >), in order 



to generate a new symmetrized basis set: {|$ M t>)} xY&uv)}- 111 (lq), Q denotes the 



g(K) elements of the group G(K) and x§ 1S the (complex) character of Q in the [K,rs\- 
representation, where K refers to one of the N allowed wave vectors in the first Brillouin zone, 
r labels the irreducible representations of the little group of K, Gl{K), and s parameterizes 
Gs- For an efficient parallel implementation of the matrix vector multiplications (MVM) 
required diagonalizing the Hamiltonian, it is extremely important, that the symmetrized 
basis can be constructed preserving the tensor product structure of the Hilbert space: 

= Nt s] rz >rs [\u) el ®\v) ph ]} . (14) 

Here u = 1, . . . , D g j K) [D 9 j K) ~ D el /g(K)], and the are normalization factors. 

19 



2. Phonon Hilbert space truncation. Since the Hilbert space associated to the phonons 



is infinite even for a finite system, we apply a truncation procedure [p9| , p7| retaining only 
basis states with at most M phonons: 

N 

{\v) P h ; rn v = ^m i)( < M}. (15) 
i=i 

The resulting Hilbert space has a total dimension D M = D 9 Jt x DK with = ^jjffi 1 , 
and a general state of the Ht JM is represented as 

l%,, s > = E E 5 -l^)- ( 16 ) 

•£t=l v=l 

It is worthwhile to point out that, switching from a real space representation to a momentum 
space description, our truncation scheme takes into account all dynamical phonon modes. 
This has to be contrasted with the frequently used single-mode approaches |34],70]. In other 



words, depending on the model parameters and the band filling, the system "decides" by 
itself how the M phonons will be distributed among the independent Einstein oscillators 
related to the N Wannier sites or, alternatively, among the N different phonon modes in 
Q-space. Hence with the same accuracy phonon dynamical effects on lattice distortions 
being quasi-localized in real space (such as polarons, Frenkel excitons,. . . ) or in momentum 
space (like polaronic superlattices, charge-density- waves,. . . ) can be studied. 

Of course, one has carefully to check for the convergence of the above truncation pro- 
cedure by calculating the ground-state energy as a function of the cut-off parameter M. 
In the numerical work convergence is assumed to be achieved if Eq is determined with a 
relative error 

. (M ) Eq(M) ~ E (M -I) a 

AE ° = eJm) - 10 • (17) 

In addition we guarantee that the phonon distribution function 

na(K) nM 

.M|2 (M)= £ £ (18) 



u=l v =l 

{m v =m} 
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which gives the different weights of the m-phonon states in the ground-state \&o), becomes 
independent of M and | C ( M )| 2 (M) < 1(T 6 . 

To illustrate the M dependences of the phonon distribution function and the ground- 
state energy, we have shown both quantities in Fig. 8 for the HtJM with e p = 4.0, Tiwq = 0.8, 
and J = 0.4. Fig. 8 proves that our truncation procedure is very well controlled even in 
the strong EP coupling regime, where multi-phonon states become increasingly important. 
Although not mentioned each time, all the results presented in this paper have been obtained 
with a sufficiently large value of M ensuring the above convergence criteria for both static 
(ground state) and dynamic (excited state) properties. 

Finally we would like to stress that it is not possible to keep all non-vanishing matrix 
elements of our ten-site Holstein t-J Hamiltonian in the memory of the present day parallel 
computers because, in spite of applying the symmetrization and Hilbert space truncation 
procedures, the total dimension of the problem is still very large. Therefore, performing the 
matrix-vector-multiply operations, we create the non-zero matrix elements of the Hamil- 
tonian simultaneously. For the example considered in Fig. 8, we have D M = 4.12 x 10 8 
(M = 15) phonons, and each MVM takes about 100 seconds cpu time on a CRAY T3E with 
126 processors. 



Algorithms for estimating spectral functions 



The numerical calculation of spectral functions, 



A°(uj) = — lim -3m 

D-l 



(iPblO* „ - - . O|^o 

u - H + E + is 



]T \(V n \Ol\%)\ 2 5[u;-(E n -E )], (19) 



n=0 



where O is the matrix representation of a certain operator O [e.g., the destruction operator 

— * 

eg a of an electron with momentum K and spin a if one wants to calculate the single hole 
spectral function (2)], involves very large sparse Hamilton matrices H. For example, for the 
EP systems studied in this paper, we have calculated spectral properties of the HtJM acting 
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in Hilbert spaces with total dimensions up to D = 10 9 . Finding all (D) eigenvectors and 
eigenstates of such huge Hamiltonian matrices is impossible, because the cpu time required 
for exact diagonalization of H scales as D 3 and memory as D 2 . Fortunately, there exist very 
accurate and well-conditioned linear scaling algorithms for a direct approximate calculation 
of A°{uj). 

In this appendix we compare the widely used spectral decoding method (SDM) ||71|| , which 
is a refinement of the standard Lanczos recursion method (LRM) ||72|| , with two more recent 
Chebyshev recursion techniques, the kernel polynomial method (KPM) |73| and the maximum 
entropy method (MEM) [[74[]. To come to the conclusion first, the MEM is found to be the 
most efficient and numerically stable method for high-energy resolution applications such 
as the lattice polaron formation problem. 

In the following we briefly review the basics of the SDM, KPM and MEM (for a more 



detailed discussion see, e.g., Refs. |35|j73| -[75| ) , and illustrate their efficiency and accuracy 
with several examples. 

1. SDM. First one calculates by iterative matrix-on-vector multiplications, 

a tridiagonal matrix [T^]^/ = (<Pi\H D \(Pii) with dimension L <ti D, starting out from \<Pi) = 
H D \<Po) — ^l^\J^ l^o); where H D denotes the Hamilton matrix acting in a .D-dimensional 
Hilbert space and \<Pq) is an arbitrary initial state having finite overlap with the true ground 
state l^o)- Applying the Lanczos recursion (|20|), the eigenvalues E n and eigenvectors \\P n ) 



of H D are approximated by those of T L |j76[, and \ty£) = Ylt=o L c niWi)i respectively. 



->L „„J |,7>L\ _ T, 



Then, having determined the ground state \&q) by the Lanczos technique, we can use 



again the recursion relation (|20|), but with the initial state 



W= , ° K) . (2D 



to determine within the SDM an approximative spectral function, 



L-l 



A°(u) = E |<o| 9 (!Pb|OtO|lPb> 6[u - {E L n - Eq)] , (22) 



n=0 
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built up by L 5-peaks. Of course, the true spectral function A°{u) has D 5-peaks. According 
to the Lanczos phenomenon the approximated spectral weights and positions of the peaks 
converge to their true values with increasing L. Some of the main problems of the SDM/LRM 
are: (i) The convergence is not uniform in the whole energy range, (ii) There exist so- 
called spurious peaks, which appear and disappear as L is increased, i.e., when the iteration 
proceeds, (iii) Without computationally expensive re-orthogonalization only a few hundred 
iterations are possible. 

2. KPM and MEM. The idea behind the Chebyshev recursion methods is to expand the 
(^-function contained in A(u) in a series of Chebyshev polynomials T m (x), 



7T\/1 — x 
with the coefficients (moments) 



A°(x) = f ] = ( + 2 f) tiT m {x)) , (23) 

m = l / 



y° = j\ixT m (x)A°(x) = (0%\T m {X)\0%) , (24) 

x = (u-b)/a,a= (E max - E min )/2(l - e), and b = (E max + E min )/2. Since the Chebyshev 
polynomials are defined on the interval [—1,1], a re-scaled Hamiltonian X = (H — b)/a has 
to be considered, whose eigenvalues x n are in the range [— (1 — e), 1 — e]. Here a small constant 
e is introduced in order to avoid convergence problems at the endpoints of the interval (a 
typical choice is e ~ 0.01 which has only 1% impact on the energy resolution |74j). The 



can be easily computed by iterative MVM. Truncating the expansion at fj,^- , the main 



difference between KPM and MEM is the way the spectral function is reconstructed from 
the remaining M moments. 

The KPM is a lin ear Chebyshev approximation to the spectrum, A (x) ~ A (x), using 
M moments only. Here the abrupt truncation of the series produces the well-known Gibbs 
oscillations which, however, can be minimized by multiplying the moments by appropriate 



damping factors W%. 



On the other hand, the MEM is a non-linear optimization procedure for calculating from 
the M known moments a representation of a (reconstructed) spectral function, A°((j)) = 
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A°(x)sm<j) with x = cos0, that corresponds to a Chebyshev series with M e ^ — F ■ M 
summands. The M e ^ moments 



p°= d(f)cos(m(f))A (4>) (25) 

•J 

are determined by maximizing the information theoretical relative entropy [[74] 



S 



(26) 



subject to the constraint that the first M moments of A°(<j)) are equal to the M known 
moments. ^4^(0) is a default model for A°((p). We choose the KPM result: A®{x) = A°(x). 
F is typically 8 ... 32. 

For the MEM (KPM) the calculation of the Chebyshev moments and not the reconstruc- 
tion of the spectral function requires the essential part of cpu time and memory. Clearly this 
is because in moment production MVM with H D is involved, whereas the reconstruction of 
the spectrum only processes the M (M) moments. 

In the numerical work we used the following algorithm producing two moments per MVM 
with H D : (i) Determine from Lanczos diagonalization the ground state |$q) and the extreme 
eigenvalues E min = E Q and E max of H D . (ii) Compute the starting conditions 

\S ) = O^\%), \S 1 ) = X\S ), (27) 

^ = (d \d ), f i° = (£ \& 1 ). (28) 

(iii) Use the recurrence relations for Chebyshev polynomials to calculate the M moments 
according to 

\S m ) = 2X|^ m _x) - \& m _ 2 ) , (29) 
f j° m = 2($ m \$ m )-n°, (30) 
f i°^ 1 = 2(£ m \£ fn - 1 )- f i°. (31) 

Note that applying this procedure only two D-dimensional vectors have to be stored. The 
algorithm is numerically very stable, so in contrast to the SDM/LRM (continued fraction 
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based technique f35fl ) thousands of Chebyshev moments can be calculated with satisfying 
accuracy. 

3. Examples. Figure 9 demonstrates the non-uniform convergence of the SDM/LRM for 
the dynamical pair correlation spectrum of the t-J model (for a discussion of the physics 
described by A a {uj) see, e.g., Refs. P5] , |75| ). Obviously, we found a much faster convergence 
in the low-energy regime. Increasing L the peak-resolution becomes better at higher energies 
as well, but now one is confronted with the re-orthogonalization problem mentioned above. 

Another typical problem concerns the proper description of gap structures. As can be 
seen from Fig. 10 (a), the gap located at the bottom of the spectrum of the d-wave pair 
correlator cannot be identified using SDM with L = 20 MVM only, because the averaged 
distance between two (SDM) £-peaks (vertical lines) is of the same order of magnitude. 
Thus the true gap structure of the spectrum remains hidden. By choosing the MEM, one 
obtains a smooth reconstruction of the whole spectrum and achieves a significant better 
resolution of the gap from the same number of MVM. This becomes particularly obvious 
by looking at the cu-integrated spectral function N a (uj). In N a (uj) the dominant low- 
energy QP peak (corresponding to the jump-like behaviour of N a (uj) at the bottom of 
the spectrum) is clearly separated from the incoherent high-energy part by an energy gap 
(with N A (u) ~ const). Of course, for both methods the energy resolution gets better with 
increasing number of moments (see Fig. 10 (b)), and consequently the SDM and MEM 
results converge on each other. However, since the dimension of the Hilbert space (D = 10 3 ) 
is rather small for the spectra shown in Fig. 10 (b), we have calculated about 10 % of the true 
peaks by using 120 MVM (on the other hand, for our EP systems having a total dimension 
of 10 9 this corresponds to only 10~ 7 of the possible excitations). Most notably, the spectral 
weight of different parts of the MEM spectrum does not change any more going from 20 to 
120 MVM. 

Since for any finite system the true spectral function will consist of a sum of delta func- 
tions, we compare the efficiency of KPM and MEM in dependence of the number of moments 
by reconstructing a fictitious reference spectrum of five 5-peaks. Figure 11 demonstrates 
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that the MEM provides a dramatic improvement in energy resolution over the KPM (typ- 
ically by a factor of four). Thus the spectral function is much better reconstructed, in 
particular when the number of MVM is small. On the other hand, an advantage of the 
linear KPM over non-linear MEM is that the KPM smoothing is uniquely determined by 
the number of Chebyshev recursions J77|] , i.e., in contrast to the MEM we have a uniform 
energy resolution oc M~ l . 

To analyze the accuracy of the KPM and MEM in more detail we have compared in 
Fig. 12 a portion of a reconstructed single hole spectrum with the true spectral function 
obtained without moment truncation, i.e., from the full diagonalization of the system. Again 
convergence to the exact data is seen to be more rapid by applying the MEM (note that 
M = AM). Moreover, the peaks in the KPM spectra are surrounded by not fully damped 
Gibbs oscillations. As can be seen from Fig. 12 (b), the MEM produces artificial peaks 
as well but they are at least one order of magnitude smaller than those for KPM. That 
means, these isolated spurious eigenvalues have extremely small spectral weight and hence 
negligible influence on the cumulative spectrum. To illustrate that these spurious peaks can 
be easily identified, we have plotted the MEM spectrum for different number of moments. 
Contrary to the true eigenvalues the spurious peaks (often labeled "ghosts" ) move about as 
M increases. 
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FIGURES 



FIG. 1. Wave-vector resolved spectral functions of one hole in the t-J (a) and Holstein 
t-J (b) models. Results are obtained for a tilted Vl~0 x y/W cluster with periodic boundary 
conditions (^-vectors are given in units of (n/5,n/5)). The integrated densities of states 
are depicted in the lower panels of (a) and (b). The units of Ag, though arbitrary, are 

— * 

nevertheless the same in all K sectors displayed in Figs. 1-3, so that the peak heights may 
be compared. 

FIG. 2. Single-hole spectral function Ag(u>) and integrated spectral weight N(u) for the 
2D HtJM at u = 0.8 for different values of e p and J. Note that the lowest peak in (c) and 
(d) contains contributions from each i?-sector. 

FIG. 3. Ag(uj) and N(uj) for the 2D HtJM in the non-adiabatic strong EP coupling 
regime (u = 3.0, e p = 4.0). 

FIG. 4. Band dispersion (a) and spectral weight factors (b) for the 2D HtJM on a 16-site 
lattice (J = 0.4). 

FIG. 5. Spectral weight factor factor Zg as a function of the EP coupling strength e p for 
the different K vectors of the ten-site lattice: (3,1) (triangles down), (2,4) (squares), (0,0) 
(circles), and (5,5) (triangles up) [in units of (tt/5, 7r/5)]. The inset shows the narrowing of 
the coherent bandwidth AE with increasing e p . 

FIG. 6. Optical conductivity in the 2D (Holstein) t-J model for J = 0.4. a re9 (oo) and 
S reg (uj) are obtained for the ten-site lattice with 15 phonons. The single-hole ground state 
has momentum K = (37r/5, 7r/5). 

FIG. 7. Optical conductivity in the 2D HtJM for J = 0.1 and c^o = 0.8. Results are given 
at e p = 1.0 (a) and e p = 4.0 in the K = (0, 0) and K = (3ir/5, 7r/5) sectors, respectively. 
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FIG. 8. Convergence of the phonon distribution function |c m | 2 (M) and ground-state 
energy E (M) (inset) as a function of the maximum number of phonons M for the 2D HtJM 
on a ten-site lattice (J = 0.4). Eq(M) is given with respect to the variational estimate of 
ground-state energy, Eq IMVLF \ obtained within the IMVLF-Lanczos approach developed 



recently by the authors [2B 



FIG. 9. Dynamical d-wave pair correlation function of the pure 2D t-J model, A a {uj) = 
En\(^ N ~ 2) \M N) )\ 2 S[co - (Ei N ~V - E ( N) )}, where A = ^ E, 5=± Qt(Q+M - PWyi)- 
Results are calculated on a 16-site lattice with two holes for J = 0.3 using SDM with 
L = 32 and 64 MVM. 

FIG. 10. SDM (thin solid lines) and MEM (thin dashed lines) results obtained for the pair 
correlation function A A (u) in the 2D t-J model with 20 (a) and 120 MVM (b) [J = 0.3, N = 
16]. Shown is also the corresponding integrated spectral weight N a (uj) = du'A (u') 
(bold lines). 

FIG. 11. Comparison of KPM (a) and MEM (b) by means of a fictitious spectrum 
consisting of five 5-peaks. The reconstructed spectral function is shown using 64 (dashed), 
256 (dotted), and 512 (solid) moments. 

FIG. 12. Depicted is a part of the excitation spectrum of a single hole with momentum 
K = (3, l)vr/5 in the 2D t-J model (D = 126). In upper panel, the KPM and MEM 
reconstructed spectral function A c (u) has been calculated using 512 and 128 moments, 
respectively. For comparison the exact peak positions and spectral weights are shown (stars). 
In the lower panel we illustrate how spurious peaks can be identified. 



32 



•a 

i 



1.0 



3 

Z 



0.0 



-5.0 



Figure 1: 



1 1 1 1 1 

(a) 

K=(5,5) 


i 




■ i i 


i 


K=(2,4) 




LL 








K=(3,l) 
















K=(0,0) 










■ i i i i i i i i 



-i 1 1 1 r~ 



0.0 



CO 



j=0.4 e =0.0 

p 



5.0 



•a 

i 



1.0 



3 

— ✓ 

Z 



0.0 



co =0.1 



-5.0 



-i 1 1 1 r~ 



0.0 



1 1 1 1 1 1 

(b) 


1. L. 


K=(2,4) 


Li / 


\k 


K=(3,l) 


l ll A ^ 




K=(0,0) 






■ i i i i i i i i I 



j=0.4 e =0.5 

p 



5.0 



to 



Figure 2: 




Figure 3: 




-8.1 
-8.3 



<N 


-8.5 - 


o 

S 


A 






-8.7 - 






w 





-8.9 
-9.1 



-9.3 



,o- 



Figure 4: 



□ □ £ p =0.2 (0..=().2 

a a e =0.5 (0 =0.8 

p o 



o... 



"O. 

-o* 



(0,0) 



(7C,0) 



K 



,D;r 

.a' 



(a) 



(it, 



0.6 



0.4 



N 



0.2 



0.0 a— 

(71,71) 



o o e p =0.0 

□ □ e p =0.2 co =0.2 

a a e p =0.5 co =0.8 

p 

/ A 



(0,0) 



(71,0) 



K 



(b) 



o 

"a\ 



a 
-a 



(it, 



Figure 5: 



T 
■ 



o.o M- 

o.o 



▲ a 

i i 

1.0 



0.8 
0.6 
3 0.4 
0.2 



— i — ■ — i — ■ — r 



0.0 



-1 , I , lJ 

().() 1.0 2.0 3.0 4.0 
£ 



co =0.8 



2.0 
8 



3.0 



4.0 



Figure 6: 




Figure 7: 




Figure 8: 




Figure 9: 




CO 



Figure 10: 




CO 



Figure 12: 




0.00 
-4.6 



-3 



0.002 



0.001 



0.000 




-4.6 



-4.4 





-4.2 
CO 



-4.0 



-3 



